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Abstract 

The mesoscopic concept is a way to deal with complex materials with an internal 
structure within continuum mechanics. It consists of extending the domain of the 
balance equations by mesoscopic variables and of introducing a local distribution 
function of these variables as a statistical element. In our case microcracks are 
modelled as penny shaped and completely characterized by their diameter and 
the unit normal to the crack surface. Two examples of crack dynamics are given as 
well as a possible definition of a damage parameter. Orientational order parameters 
(fabric-alignment tensors) are defined and balance like dynamic equations for them 
are derived. 



1. A MODEL OF MICROCRACKS 

Macroscopic failure of brittle materials is initiated by the propagation of micro- 
cracks. In a simple model the microcrack is described as a flat, rotation symmetric 
surface, a so called penny shaped crack. In addition we make here the following 
simplifying assumptions: 

(1) The diameter of the cracks is much smaller than the linear dimension of the 
continuum element. Under this assumption the cracks can be treated as an 
internal structure of the continuum element. The cracks are assumed small 
enough that there is a whole distribution of crack sizes and orientations in 
the volume element. 

(2) The cracks are fixed to the material. Therefore their motion is coupled to 
the motion of representative volume elements. 

(3) The cracks cannot rotate independently of the material, i. e. if they have 
a nonzero rotation velocity at all, this rotation velocity is determined by 
the antisymmetric part of the time derivative of the deformation gradient 
of the surrounding material and it does not depend on crack length and 
orientation. All cracks within a volume element move and rotate with the 
same velocity. 

(4) The number of cracks is fixed, there is no production of cracks, but very 
short cracks are preexisting in the virgin material. 

(5) The cracks cannot decrease area, but can only enlarge, meaning that cracks 
cannot heal. 

To summarize our model the microcrack is characterized by a unit vector n repre- 
senting the orientation of the surface normal and by the radius I of the spherical 
crack surface. These parameters will be taken as the additional variables in the 
mesoscopic theory. 
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2. Different approaches to damage mechanics and the mesoscopic 

CONCEPT 



There are two principally different possibilities to deal with complex materials 
within continuum mechanics: The first way is to introduce additional fields de- 
pending on position and time. These fields can be any kind of internal variables 
|l| 0], or damage parameters || |J, and damage tensors (fabric tensors) || H{. In 
damage mechanics such additional macroscopic variables have been introduced in 
many different cases of materials with internal structure like liquid crystals j?], el, 
polymer solutions || [l(]] and others. The other way is a so called mesoscopic theory. 
The idea is to enlarge the domain of the field quantities by an additional variable, 
characterizing the internal degree of freedom connected to the internal structure 
of the material. Field quantities are introduced, which are defined on an enlarged 
space IjXlfXM. The manifold M is given by the set of values the internal degree 
of freedom can take. In our case the internal degrees of freedom are the different 
sizes / and orientations n of microcracks. We call this way of dealing with the in- 
ternal structure of complex materials a mesoscopic theory, because it includes more 
information than a macroscopic theory on R^ x R t , but less than a microscopic one 
on the molecular level. The domain of the mesoscopic field quantities x R t x M 
is called mesoscopic space. 

Macroscopic quantities are calculated from mesoscopic ones as averages over 
crack sizes and crack orientations. The spatial distribution of cracks is not relevant 
in the sense that the resulting macroscopic quantities are still field quantities de- 
pending on position and time. For a treatment of the spatial distribution of cracks 
and a possible coarsening process see Q . 

In contrast to spatial averaging introduced in [l^] , a nonlocal generalization of 
the classical Weibull-theory, the averages in mesoscopic theory are local in space. 
They are averages over different microcrack sizes and orientations in a volume 
element. 

Scaling properties |13 are completely out of scope of the whole mesoscopic the- 
ory. They result from microscopic statistical considerations. Statistical theories 
of fracture describe the breakdown of material as a second order phase transition 
14, 13 , as well as a first order phase transition |H| [l6|, [l7j. 

We will apply now the mesoscopic concept to damaged material with micro- 
cracks. The crack length can take values between a minimal length l m of the 
smallest preexisting cracks and a maximal length Im, which is limited by the linear 
dimension of the continuum clement. The orientation of the unit vector n can be 
given by an element of the unit sphere S 2 . Therefore in the example of microcracks 
the manifold M is given by [l m , Im] x S 2 . The change velocities of the mesoscopic 
variables I and u .— fi are defined in such a way that for At — > we have 



(2.1) l(t + At) = l(t) + lAt , n(t + At) = n(t) + uAt 

at later times t + At. The rotation velocity u and the length change velocity I are 
the components in spherical coordinates of the crack velocity vi introduced in Jl8| . 
In this previous paper jl^] the set of additional mesoscopic variables n and I was 
called directional variable. 

Beyond the use of additional variables the mesoscopic concept introduces a sta- 
tistical element, the so-called mesoscopic distribution function. In our case this is 

2 



a distribution of crack lengths and orientations in the continuum element at posi- 
tion x and time t, called here crack distribution function (CDF). The distribution 
function is the probability density of finding a crack of length I and orientation n 
in the continuum element. 



3. Mesoscopic balance equations 

Now fields as mass density, momentum density, angular momentum density, and 
energy density are defined on the mesoscopic space. For distinguishing these fields 
from the macroscopic ones we add the word "mesoscopic". In addition we introduce 
the crack number density N as an extensive quantity. The mesoscopic crack number 
density N(l, n, x, t) is the number density, counting only cracks of length I and 
orientation n. For this crack number density there is a balance equation too, as it is 
an extensive quantity. The crack number density can be prescribed independently 
of the mass density, although the motion of cracks is coupled to the motion of 
surrounding material in our model. Therefore we distinguish here between the two 
fields mass density p and number density TV, although they have the same equation 
of motion and where not distinguished in an earlier paper Jl8| . 

3.1. Definition of the distribution function. Due to its definition as probabil- 
ity density the distribution function is the number fraction 

(3 . 2) /(i ,„„ t) .toi , 

in volume elements, where the number density N(x,t) is non-zero. Here N(x,t) 
is the macroscopic number densit y of cracks of any length and orientation. As the 
distribution function in equation (|3.2[ ) is not well defined if N(x, t) = 0, we define 
in addition that in this case f(l, n, x, t) = 0. As there is no creation of cracks in our 
model the distribution function will be zero for all times in these volume elements. 
In all other volume elements with a nonzero crack number it is normalized 

(3.3) / / f(l 7 n,x,t)l 2 d 2 ndl = l . 

Jim JS 2 

3.2. Balance equations of mass, momentum, angular momentum, and en- 
ergy. For the mesoscopic densities local balance equations have been derived from 
the macroscopic global ones |E| ^0|, [2ll [l^]. The macroscopic balance equations 
express the fact that the extensive macroscopic quantities within a region G can 
change due to a flux over the boundary dG and due to production and supply 
within G. This results in the general form of a global balance equation 

(3.4) f Xd 3 xdnl 2 dl = [ 4> x (-)da + [ Z x (-)d 3 xdnl 2 dl . 

dt Jg JdG Jg 

A generalized Reynolds transport theorem on the mesoscopic space, analogous to 
the one in p2| , is used to transform the time derivative, and a generalized Gauss 
theorem is applied. The boundary dG of G consists of parts in position space, 
in orientation space, and in the length interval. In regular points of the contin- 
uum we get the general form of a local mesoscopic balance equation JlS| with the 
abbreviation (•) = (I, n, x,t): 

^X(-)+V x -b(-)X(.)-5(.)]+V n -[ W (-)X(.)-i2(-)]+i| (l%)X(-) Ri(-)) = £(• 



where R and Ri are the non-convective fluxes over the orientational and length part 
of the boundary of G, and G is a region in M 3 x S 2 x [l m , Im\- The derivative with 
respect to the mesoscopic variable {I, n) is represented in spherical coordinates. In 
the derivation of the local balance equation it has been supposed that there is no 
flux over the boundary of the total mesoscopic space: 



(3.5) 



,.00 * A M 

/ / / V-cpx(-)d 3 xdnl 2 dl = 

J-oo Js 2 JU 



Otherwise such a non-zero flux term (3.5) could be interpreted as an additional 
source term on the right hand side of the equation. 

Explicitly we have: 
Balance of mass 



(3.6) ^g(-) + V x -{g(-)v(x,t)} + V r , 



0. 



Balance of momentum 
d_ 
dt 

u(x, t)g(-)v(x, t) 



+V„ 



[g(-)v(x,t)]+V a 



T T (.) 



[v(x,t)g(-)v(x,t) 
1 d 



I 2 dl 



l z lg{-)v{x,t) -r(-) 



(3.7) 



= £?(•)/(•)• 

Here /(•) is the external acceleration density, t T (■) the transposed mesoscopic 
stress tensor, and T T (■) the transposed stress tensor on orientation space (non- 
convective momentum flux in orientation space), r(-) is the momentum flux vector 
with respect to the crack length variable. Wc introduced already the assumption 
that the material velocity v and the rotational velocity u are the same for cracks 
of all orientations and lengths. 



Angular Momentum 

The balance of angular momentum has to be taken into account as an additional 
equation independent from the balance of momentum, because there is an internal 
angular momentum due to crack rotations. The total angular momentum 

(3.8) S(x,t) := x x v(x,t) + s(-), 

is the sum of the moment of momentum and the internal angular momentum. 

|[f?(W)] + V. • [v(x, t)g(-)S(-) (x x T(-)) T - H T (.) 

+V„ ■ [u(x,t)g(-)S(-) (x x r(-)) T - tv t (■)] +~ (l 2 ig(-)S(-) «(•) 

(3.9) =g(-)x x fc(-) + g(-)g. 

Here n x g is the vector of couple forces (acting on crack orientation), the tensor 
II is the surface torque, and 7r is the analogue with respect to orientation, and u> 
is the analogue with respect to crack length. These constitutive quantities appear 
in the non-convective fluxes in position space, orientation space, and in the length 
interval, respectively. This equation can be simplified with the assumptions that 
the material velocity and the rotation velocity depend only on position and time 
v(x, t) and u(x, t). 
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However, the spin is only relevant, if the model allows for crack rotations inde- 
pendently from the rotations of material elements, and this is not the case in our 
simplified example dynamics. 



Similarly the balance of energy can be given, which is omitted here and can be 
found in fL8f . In all balance equations in addition to the flux with respect to the 
position variable there appear additional flux terms with respect to the additional 
mesoscopic variables crack orientation and length. 

Balance of crack number 
In our model the cracks move together with the material clement. Therefore their 
flux is the convective flux, having a part in position space, a part in orientation 
space, and a part in the length interval. There is no production and no supply of 
crack number. Therefore we have for the crack number density N: 



(3.10) ^N(-) + V, • {#(>(*, <)} + V„ • {N(-)u(x, t)} + ~ (l 2 iN(-)) 



= 0. 



In a fixed volume element this crack number density is proportional to the mass 
density, and therefore these two fields where not distinguished in an earlier paper 



18|. 



We obtain from the mesoscopic balance of crack number density (3.1C) a balance 
of the CDF f(l,n,x,i), by inserting its definition (|3T 



d 

— f(l, n, x, t) + \7 x - (v(x, t)f(l, n, x, t)) - 

I Q , . 

V„ • (u(x,t)f(l,n,x,t)) + ^ — (l 2 lf(l,n,x,t) 
-f(l,n,x,t) ( d 



N(x,t) Kdt +v{x,t).V,)N(x,t) 

_ -f(l,n,x,t)dN(x,t) _ n 
~ N(x,t) dt ~ U ' 

The right hand side is equal to zero, as for the co-moving observer the total num- 
ber of cracks in a volume element does not change in time, as in our model no 
cracks are created. In our model all cracks in a volume element move with the 
translational velocity of the volume element v(x,t) and rotate with the velocity 
u = V x v(x,t). Therefore the first three terms on the left hand side can be 
summarized as a co-moving time derivative of the distribution function (the time 
derivative of an observer moving with the material elements) with the abbreviation 

J/ ' 

d 

— f(l, n, x, t) + v(x, t) ■ V x f(l, n, x, t) + 
(3.12) u(x,t).V n f(l > n,x,t) = dCfil '£ X > t) . 

If we assume in addition an incompressible motion W x ■ v = 0, we end up with the 
equation of motion for the CDF: 

( , 13) ^M + i|(, 1/(1 ,„,,, i)H 



This is not yet a closed differential equation for the CDF as long as no expression 
for the length change velocity of the crack I is given. An example of such a closed 
equation will be discussed later. 

Macroscopic quantities are obtained from mesoscopic ones as averages with the 
CDF as probability density: 

plM r 

(3.14) A(x,t)= / / A(l,n,x,t)f(l 7 n 7 x,t)d 2 nl 2 dl 

Jim Js 2 

Entropy balance 

Besides these mesoscopic balances the entropy balance is necessary for introducing 
the second law of thermodynamics. Because the production of mesoscopic entropy 
is not necessarily positive for each crack length and orientation, the entropy balance 
is only interesting in its macroscopic form. 

d 

— [g(x, t)T}(x, t)\ + V x ■ [g(x, t)r](x, t)v{x, t) + <f>(x, t)\ = 

(3.15) = g(x,t)a(x,t) 

(ji = specific entropy density <p — entropy flux density a = entropy production 
density). The second law is expressed by the dissipation inequality 

(3.16) <r(x,t)>0 . 

The set of balance equations is not a closed system of equations, constitutive 
equations for mesoscopic quantities are needed. The domain of the constitutive 
mappings is the state space; here a mesoscopic one. There are the possibilities that 
the mesoscopic state space includes only mesoscopic quantities, or that it includes 
mesoscopic and macroscopic quantities, and there are examples where such mixed 
state spaces cannot be avoided Q. (For instance in the case of liquid crystals 
the macroscopic alignment tensor is included in a mesoscopic state space. This is 
necessary to account for the orienting mean field of surrounding ordered particles. 
Otherwise it is not possible to describe the phase transition from the isotropic phase 
to the ordered liquid crystalline phase.) Constitutive equations have to be such that 
the second law of thermodynamics is fulfilled by any solution of the macroscopic 
balance equations with the constitutive equations inserted p3|. This requirement 
restricts the possible constitutive functions. 

Finally, even for the exploitation of the dissipation inequality, which is possible 
only on the macroscopic level, the choice of variables can be motivated by the 
mesoscopic background |^5|, [2(| . A relevance of these variables could not be guessed 
from a purely macroscopic theory. 



3.3. Damage parameter and order parameters. 



3.4. Definition of a damage parameter. The damage parameter is introduced 
as a macroscopic quantity growing with progressive damage in such a way that 
it should be possible to relate the change of material properties to the growth of 
the damage parameter. We define the damage parameter as the fraction of cracks, 
which have reached a certain length L. The idea is that cracks of this and larger sizes 
considerably decrease the strength of the material, and therefore their fraction is a 
measure of the damage. This idea is related to the slender bar model of Krajcinovic 
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|5|, where the damage parameter is introduced as the number of 'broken bars' in 
the sample . 



(3.17) D{x,t)= [ [ f{l,n,x,t)d 2 nl 2 dl . 

Jl Js 2 

In this definition of the damage parameter the possibility of cracks of any length 
(Im ~~ * 00 ) is included. This is consistent with many possible laws of crack growth, 
where the crack does not stop growing. 

More sophisticated definitions, taking into account the orientational distribution 
too, are possible and will be discussed elsewhere. Another measure of damage, 
which could be introduced is the average crack length (27jj. 

3.5. Length order parameters. From the mesoscopic distribution function two 
different kinds of moment series can be built because of the dependence on crack 
length and on crack orientation: We can introduce moments of the distribution 
function with respect to crack length: 

/•Im 

(3.18) ^ f{l,n,x,t)P k (l)l 2 dl =: p k (n,x,t) 



where P k (l) are polynomials being orthogonal with respect to the measure l 2 dl: 

Am 



(3.19) / P % (l)P 3 {l)l 2 dl = Sij 



The moments introduced in equation (3.18) still depend on crack orientation. 
Averaging over all orientations gives macroscopic fields, the length order parame- 
ters: 

(3.20) 7r*(x,t)= / P k (n,x,t)d 2 n . 

Js 2 

In the following we will investigate the moments of the distribution function with 
respect to crack orientation. 

3.6. Orientational order parameters. We can introduce the following set of 
alignment-fabric tensors of successive tensorial order 

Am r 

(3.21) a {k) (x,t) := / f(l,n,x,t) no... on l 2 dld 2 n 

Jl m Js 2 ' - ' 

where tt? denotes the symmetric irreducible part of a tensor ]2Sj] . Remarkable, 
that only the even order tensors appear in the series because the microcracks arc 
represented by axial vectors, the unit normal to the crack surface, i.e. n and — n are 
not distinguished. Due to this symmetry all odd order moments vanish. The tensors 
defined above are macroscopic quantities. We want to call them alignment-fabric- 
tensors. Originally tensorial damage parameters where introduced on a purely 
statistical ground, without a mesoscopic foundation and where called 'fabric tensors 
of the second kind' in damage mechanics (see Kanatani Q or Krajcinovic [El). The 
alignment- fabric-tensors represent the orientational distribution of microcracks, but 
do not take into account their lengths. They have to be distinguished from the scalar 
damage parameter which is a measure of the growth of cracks. These alignment- 
fabric tensors form a whole set of internal variables in the sense of thermodynamics. 
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The alignment-fabric tensors are a measure of the deviation of the crack orien- 
tation distribution from isotropy. They are all zero, if all crack orientations are 
equally probable, and at least some alignment-fabric tensors are nonzero in case 
of anisotropic distributions. The orientation distribution of cracks and therefore 
the alignment-fabric tensors become important in the dynamics of the crack distri- 
bution. There are usually the specimen geometry and loading conditions rotation 
symmetric around an axis d (uniaxial conditions). It is reasonable to assume that 
also the distribution of crack orientations is rotation symmetric around the same 
axis d. Then, for symmetry reasons, all alignment-fabric tensors of different orders 
can be expressed in terms of scalar orientational order parameters and the unit 
vector d in the following way: 

(3.22) a (fc) - g(fc) o . . . Q J (A; = 2,4, ...) , 

fe 

where the order parameters are one in case of total alignment (the microcracks 
stand parallel) and zero for randomly oriented cracks. 

3.7. Equations of motion for the alignment-fabric tensors and for the 
damage parameter. In general a coupled set of equations of motion for the 
alignment- fabric tensors of different order can be derived from the differential equa- 
tion for the crack distribution function by taking moments of this equation, i.e. 
multiplying with the dyadic product n o . „ o n 1 and integrating over all orienta- 

k 

tions n S S 2 . This set of equations is analogous to the differential equations for 
the alignment tensors in liquid crystal theory Jl9| and will be discussed elsewhere 
in more detail. In general the equations for the different tensor orders are coupled. 

As in our model all cracks in a volume element have the same angular velocity, 
namely that of the surrounding material, this set of equations simplifies to a set of 
very special balance type equations without production, and without non-convective 
flux, which are not coupled: 

f) n (k) 1 1 

(3.23) — \-v(x, t) • Va (fc) + - (V X v) ■ a (fe) - -d k > • (V X v) = 

(3.24) or a (fe) = 

for any tensor order k. This special form arises due to the model assumption of a 
fixed crack number and in addition cracks not moving and rotating independently 
of the surrounding material. Therefore for an observer co-moving with the material 
the orientation distribution and the alignment-fabric tensors do not change. These 
equations are the equations of motion for the internal variables, which have to be 
postulated in a purely macroscopic theory. For our simplified crack dynamics the 
dynamics of the alignment-fabric tensors is not independent of the motion of the 
material elements. Therefore the change of the alignment tensors in time is not 
relevant to be considered in our simplified model, as it is completely determined 
by the motion of the surrounding material. However, the situation is different for 
other, more complicated crack dynamics. In any case, even if the dynamics of the 
tensorial damage parameters is not interesting, the orientation distribution itself is 
relevant, because of the dependence of the effective stress on crack orientation (see 
below). This effective stress determines the dynamics, as it appears for instance in 
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the Griffith-criterion for the onset of growth, and it also appears in the expressions 
for the length change velocity discussed in the examples below. 
Orientation dependence of the effective stress 



In an experiment with uniaxial tension a applied to the sample (see figure 4.2) the 
stress component cr n , normal to the crack surface, depends on crack orientation. 

Let us assume that in an experiment a uniaxial tension a is applied along the 
z-direction. Then the stress component in the direction n, acting on a crack surface 
with surface unit normal n is 

(3.25) cr e ff = a(e z ■ n) 2 

where e z is the unit vector in z-direction. This dependence of the effective stress 
on crack orientation leads after averaging over all orientations to 

/ a e fff(l,n,x,t)d 2 n= / a(e z ■ n) 2 f(l,n,x,t)d 2 n 
Js 2 Js 2 

[ a ( nn — —8) f(l,n,x,t)d 2 n + -8 j af(l,n,x,t)d 2 n 
Js 2 V 3 / 3 Js 2 



e z e z 



(3.26) 



3 



e z e z 



where a zz is the zz-component of the second order alignment-fabric-tensor a. This 
dependence of the effective stress on the alignment-fabric-tensor leads to a depen- 
dence of the crack dynamics (for instance the critical length) on the orientational 
order. Thus macroscopic equations of motion of damage parameters depend on the 
orientational order characterized macroscopically by the alignment-fabric tensors. 
Hence it would be interesting to study the dynamics of the alignment-fabric-tensors, 
too. 

3.8. Differential equation for the damage parameter. Differentiating the def- 
inition of the damage parameter equation (3.17) with respect to time we get the 
following differential equation for the damage parameter: 

dD(x,t) _ d r h 
~dt] L 



dt 



f(l, n, x, t)d 2 nl 2 dl 



— /(/, n, x, t)l z + f(l, n, x, t)2ll d 2 ndl 
52 v dt 

d_ 

di 



(l 2 f(l,n,x,t)lj + f(l,n,x,t)2li) d 2 ndl 



(3.27) 



l 2 f(l,n,x,t)l +2 / f(l,n,x,t)lld 2 ndl 
J L Jl Js 2 

The differential equation for the damage parameter depends on the crack dis- 
tribution function itself, and therefore also on the initial crack distribution, and it 
also depends on the dynamical equation for the crack length. 

4. Examples of closed differential equations for the distribution 

function 

Some model on the growth velocity of a single crack is needed in order to make a 
closed differential equation for the length and orientation distribution function out 
of equation ( 3.1 1| ) . Two different dynamics of crack extension from the literature 
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will be given here as examples. In the second example we suppose that for a given 
load not all cracks start growing, but only cracks exceeding a certain critical length 
l c , which is given by the Griffith-criterion. As in many examples of a crack length 
change dynamics the cracks do not stop growing, but extend infinitely, in all these 
cases the maximal crack length has to be set to Im = oo. However, when the cracks 
become macroscopic their growth dynamics becomes more complicated (showing 
for instance branching) than our example dynamics. 

4.1. Mott extension of Griffiths energy criterion including a kinetic crack 
energy. When the cracks are growing the system has a kinetic energy due to the 
growth by virtue of the inertia of the material surrounding the separating crack 
surfaces. This extension of the original Griffith energy concept (see below) by a 
kinetic energy term goes back to Mott j29): A kinetic energy term is added to 
the sum of the crack surface energy and the elastic deformation energy of the 
surrounding elastic material, and the crack length is such that the total energy of 
the system is constant. Two different loading conditions are especially interesting: 
fixed loading ('dead weight') and 'fixed grips' conditions. In both experiments 
uniaxial symmetry is assumed. In the first case a constant force is applied to 
the ends of the specimen, leading to a tensile stress. In the second case a fixed 
displacement is prescribed at the outer boundaries of the specimen. For these two 
loading conditions, requiring a constant total energy and an argument based on 
geometrical similarity, the following expressions for the crack length change velocity 
have been derived ([|30[p93): 
'Dead weight': 

(4.28) 1 = l T (l-j 

where It is the so called terminal velocity, not depending on crack length, but on 
the applied load (Jeff, and therefore on crack orientation, is the initial crack 
length. 
'Fixed grips': 

(4.29) i = i 

where the parameter a is defined as 

(4.30) 8 ^ 




A 

It is the ratio of the initial crack area to the surface area A of a cross section of the 
specimen. In the 'fixed grips' geometry the crack extension might stop again after 
a certain growth. This can be understood, because of the increase in compliance 
associated with crack extension in a finite specimen. This leads to a diminishing 
applied force and decreasing tendency of the crack growth. 



From the mesoscopic point of view the growth laws equations (4.30) as well as 



(f.28j) are mesoscopic constitutive equations relating the length growth velocity I 



to the external load in a material dependent manner. 

In both loading conditions discussed here the crack velocities have been derived 
for single cracks. If we apply these growth velocities in our differential equation for 
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the length distribution function, eq. (3.11) this means that we neglect interaction 
between cracks. However, crack interactions can be taken into account by more 
sophisticated expressions for the length change velocity. 

Inserting the length change velocities of the previous section into the differential 
equation for the crack distribution function, and integrating over all orientations 
leads to the following closed differential equations: 
'Dead weight': 



(4.31) 



df(l,x,t) 
dt 



The parameter It depends on the effective load cr e ff and therefore on the second 
order alignment-fabric tensor. 
'Fixed grips': 



df(l,x,t) 



(4.32) 



1__5 



( 



l 2 f(l,x,t)h 



V 



1 




4.2. Griffith-criterion for the onset of growth. The criterion for the cracks to 
start growing adopted in the example is the energy criterion introduced originally 
by Griffith ||. According to Griffith |l| there is a criticality condition for the 
crack growth to start, and for cracks larger than a critical length there is a velocity 
of crack growth L From energetic considerations Griffith (3^] derived a critical 
length of cracks with cracks exceeding this length starting to grow. This critical 
length is given by : 



(4.33) l c = — , 

where if is a material constant, and a n is the stress applied perpendicular to the 
crack surface. It is assumed that a stress component within the crack plane does 
not cause crack growth. For cracks smaller than the critical length l c the energy 
necessary to create the crack surface exceeds the energy gain due to release of 
stresses. 

4.3. Rice-Griffith dynamics. A possible crack dynamics, taking into account the 
criticality condition of Griffith is derived from a generalization of the Griffith en- 
ergy criterion on thermodynamic grounds, introducing Gibbs-potential |32]j , which 
includes the stress normal to the crack surface and crack length as variables. The 
resulting crack evolution law has the form 

(4.34) i=-a + f3a 2 l for I > l c , 

(4.35) 1 = for I < l c 

with material coefficients a, and f3. In case of a constant time rate of the applied 
stress, er = v ff t, it results: 

li 



(4.36) l = -a + (3v 2 Jt 2 for I > l c , 

(4.37) 1 = for I < l c . 

v a is the time derivative of the applied stress normal to the crack surface. The 
dependence of this normal stress on crack orientation leads to the following orien- 
tation dependence of the dynamics: 



(4.38) I = -a + f3v 2 a lt 2 (e z ■ n) 4 for I > l c , 

(4.39) 1 = for I < l c , 

where v a o is the change velocity of the stress applied in the z-direction. 

After averaging over all orientations this orientation dependence leads to a de- 
pendence on the fourth moment J„ 2 nnnn fcPn of the distribution function. 

This dynamics also includes a criticality condition for the crack to start growing. 

With this model for the length change velocity we end up with the following 
differential equation for the distribution function: 

df(l,n,x,t) 1 d , 2 . („\*iJr\\ 

(4.40) for I > l c 

df(l,n,x,t) 

4.41) M ' ' ' =0 for I < l c . 

Solutions of this differential equation will are discussed in M. 



5. Conclusions 

In the mesoscopic description we have introduced mesoscopic fields, defined on an 
enlarged space including crack size and orientation. Averages over crack sizes and 
orientations, i.e. macroscopic quantities are calculated with a distribution function 
/. The differential equation for this distribution function was derived from the 
mesoscopic balance equations and crack growth law for the single crack. Different 
such crack growth laws from the literature were discussed. 

Macroscopic quantities accounting for the progressive damage have been defined 
as integrals calculated with the distribution function. These are scalar damage 
parameters, like for instance the average crack length, and Fabric-alignment tensors. 
For these different scalar and tensorial damage parameters equations of motion have 
been derived. The time evolution of Fabric- alignment tensors will be of special 
importance under biaxial loading conditions. 

The equations of motion for the damage parameters can be compared to the 
evolution equation in phase field models (or in Landau-thoery of phase transitions. 
In phase field models an additional wanted field, the phase field is introduced. 
The form of the equation of motion, often in the form of a conservation law is 
postulated |3^, This phase field can be compared to the damage parameter 
introduced here, and in the non-unilateral case also to the Fabric-alignment-tensor. 
The equation of motion for the damage parameter is of the same type. It is a 
special form of a balance equation, here with a zero flux term, because spatial 
inhomogeneities were not taken into account. However, this form of equation of 
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motion has not been postulated here, but derived from mesoscopic considerations, 
i.e. mesoscopic balance equations. 
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